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ABSTRACT 

We present the 250, 350, and 500 /im detection of bright submiUimeter emission in the direction of 
the Bullet Cluster measured by the Balloon-borne Large- Aperture SubmiUimeter Telescope (BLAST). 
The 500 fim centroid is coincident with an AzTEC 1.1 mm point-source detection at a position close to 
the peak lensing magnification produced by the cluster. However, the 250 /im and 350 fim centroids are 
elongated and shifted toward the south with a differential shift between bands that cannot be explained 
by pointing uncertainties. We therefore conclude that the BLAST detection is likely contaminated 
by emission from foreground galaxies associated with the Bullet Cluster. The submiUimeter redshift 
estimate based on 250-1100 /im photometry at the position of the AzTEC source is Zphot = 2.9^q'3, 
consistent with the infrared color redshift estimation of the most likely IRAC counterpart. These flux 
densities indicate an apparent far-infrared luminosity of Lfir = 2 x 10^'^ Lq. When the ampliflcation 
due to the gravitational lensing of the cluster is removed, the intrinsic far-infrared luminosity of the 
source is found to be LpiR 10^^ Lq, consistent with typical luminous infrared galaxies. 
Subject headings: submiUimeter — galaxy clusters 



1. INTRODUCTION 

A decade of advances in submiUimeter surveys has 
led to a broad understanding of the galaxies responsi- 
ble for the high-re dshift component of the cosrnic in- 



frared back ground (i Smail et al. 



1997t iBarger et al 



1998 



2004 



2007 



Hughes et al. 1998; W ebb et al.ll2003l: iGreve et al 
Ppve et al. 2006; Coppin et al."2006VBert oldi et al.l 
Knudscn et al. 2008; Devlin et al. 2009). These massive 
submiUimeter galaxies (SMGs) are considered to be the 
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products of galaxy mergers and are associated with a vi- 
tal stage of galaxy evolution. With typical far-infrared 
(FIR) luminosities > 10^^ L0, they are presumed to 
be the high-redshift counterparts to (ultra) luminous in- 
frared galaxies (LIRGs, ULIRGs). The source of this en- 
ergy is generally attributed to powerful starbursts, sug- 
gesting star formation rates of 100-1000 MQyr~^. Ap- 
proximate ly half of these galaxi e s are located at 1.9 < 
z < 2.9 (jChapman et al.l l2005t [Aretxaga et"all 120071) ^ 
constituting a significa nt component of the star fo rma- 
tion rate at this epoch (jPerez-Gonzalez et al.l[2005l ). 

The increased mapping speeds of current submiUime- 
ter instruments have facilitated rec e nt large-area blank- 
field surveys (iSiringo et all [200l iDevlin et all [20091 : 
lAustermann et al.ll2009i r These surveys provide a large 
statistical sample of bright SMGs. However, resolu- 
tion restricts the depth to the confusion limit and they 
are unable to probe the abundant faint end of the 
number counts. Observations directed toward massive 
galaxy clusters take advantage of gravitational lensing 
to penetrate the confusion limit an d sample intrinsi- 
cally fainter and higher-redshift SMGs (iSmail et al.|[l997l : 
iZemcov et £01120071: iKnudsen et al.ll2008f K 

An extremely bright millimeter source was recently re- 
ported from a survey of the Bullet Cluster conducted 
with the 1.1mm AzTEC receiver on the Atacama Sub- 
miUim eter Telescope Experiment (ASTE) (| Wilson et al.l 
120081) . Its millimeter properties are given in Table [TJ 
The high apparent flux and close proximity to a crit- 
ical line of magniflcation on the lensing model moti- 
vated the authors to conclude that the source is likely 
to be a background galaxy amplified by the foreground 
mass. This AzTEC source (MMJ065837-5557.0) is co- 
incident with a LIRG detected by IRAC (which possibly 
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has multiple infrared images of the galaxy), with a red- 
shift z^^^t = 2.7±0.2 (IWilson et al.llfinlirTonzalez et all 
l2009h . 

We present the 250, 350, and 500 /im measurements 
of a bright source (BLAST J065837-555708) located 
in the direction of the Bullet Cluster. With signifi- 
cances 6-(T, 5-cr, and 9-a at 250, 350, and 500 /im, re- 
spectively, this detection is the highest signal-to-noise 
source in the 0.25 deg^ map of the region taken with the 
Balloon-borne Lar ge-Aperture SubmiU imeter Telescope 
(BLAST) in 2006 (jPascale et al.ll2008[ ). Details of the 
observations and data reduction techniques are presented 
in Section [21 Extracted photometry as well as a photo- 
metric redshift estimate and constraints on the intrinsic 
properties of the millimeter source are given in Section[3l 
Finally, our conclusions are summarized in Section 3] 

2. OBSERVATIONS AND REDUCTION 

BLAST conducted an 11 day high-altitude flight from 
McMurdo Station, Antarctica, in December of 2006 
(BLAST06). The telescope has a 1.8 m primary mirror 
and three large-format bolometer arrays that image the 
sky simultaneously at 250, 350, and 500 /im over a 14' x 7' 
field of view. The telescope point spread function (PSF) 
at each wavelength is best fit by Gaussians with full- 
widths at half-maximum (FWHM) 36", 42", and 60" at 
250, 350, and 500 /im, respectively. BLAST scans the sky 
quickly (~ 0.1 deg s^"'^) in azimuth, while drifting slowly 
in elevation. Details o f the BLAST instrument can be 
found in Pascale et all 12008). The BLAST wavelengths 
were selected to sample the peak of the spectral energy 
distribution for high-redshift starburst galaxies as well as 
cooler star-forming regions within our own Galaxy. 

A large portion of the BLAST06 flight (> 150/220 
hours) was dedicated to conducting large-area, wide and 
deep blank-field surveys to identify and characterize the 
high-redshift population of starburst galax ies that com- 
prise the FIR extra-galactic background (jDevlin et al.l 
I2009f) . Two smaller fields were also targeted which are 
known to contain massive galaxy clusters, including a 
1.1 deg'^ map around Abell 3112 and a 0.25 deg'^ map 
around the Bullet Cluster. The Bullet Cluster field was 
visited three times during the flight, with a total observ- 
ing time of 2.0 hours. 

The raw BLAST06 time-stream data are redu c ed us - 
ing a common pipeline detailed in ' Pascale et al.l (|2008l ) 
and lTruch et a l. (2008). The absolute calibration of the 
BLAST06 data is determined based on observations of 
the star VY CMa with 10%-13 % uncertainty wh ich is 
highly correlated between bands (|Truch et al.ll2009f ). Sig- 
nal and variance maps of the Bullet Cluster with 10" pix- 
els are generated using the mapmaking algorithm Opt- 
Bin (Pascale et al., in prep). The rotation of the sky 
between visits to the field provided '^10° cross-linking 
between azimuthal scans. This shallow angle is insuffi- 
cient to constrain low frequency drifts in the detectors 
and preliminary maps suffered from significant striping 
along the telescope scan direction. In order to suppress 
these artifacts, the time-stream data are filtered by the 
map maker, whitening the map on scales > 16'. The 
appropriate PSF for each band is produced by stack- 
ing maps of VY Cma, co-added with the the proper sky 
rotation, and filtered in the same manner as the final 
maps to enable accurate source extraction and calibra- 



tion^. A catalog of sources for each band is compiled 
using a source-finding algorithm which selects the peaks 
in a smoothed map produced by the noise-weighted con- 
volution of the image with the PSF. The peak positions 
in each band associated with BLAST J065837-555708 
are given in Table [2] 

The telescope pointing solution is reconstructed post- 
fiight by combining 100 Hz rate gyroscope data with ~ 
1 Hz star camera solutions w hich provide an absolute 
reference (jPascale et al.l l2008| ) . The star cameras are 
mounted to the same rigid structure as the telescope. 
Nevertheless, tilting the structure in elevation results in 
small changes to the angular offset between the bore- 
sights. To correct this effect, bright compact objects 
were targeted as pointing calibrators at different eleva- 
tions throughout the fiight. The resulting 1-cr pointing 
uncertainty over the elevation range of observed calibra- 
tors (30°-50°) is - 5". Observations of the Bullet Clus- 
ter were performed at high elevation (^ 54°), outside 
the range of our pointing calibrators. They were also 
conducted early in the fiight, before the telescope focus 
had been finalized. While these factors lead us to suspect 
larger absolute pointing uncertainties for sources in the 
Bullet Cluster field, the optical design of BLAST does 
not allow for relative offsets between the bands. 

The absolute pointing offset for the BLAST maps of 
the Bullet Cluster is determined by stacking them against 
positions of 24 /im sources. This 24 /xm catalog con- 
tains sources with a significance > 3-(j extracted from 
the Spitzer MIPS observations of the field usi ng the 
APEX source-finding algorithm (|Makovoz et all [2002). 
At 250 /im, this analysis produces the highest signal-to- 
noise, indicating a 27" astrometry correction required 
for the BLAST maps. The corrected peak positions of 
BLAST J065837-555708 are given in Table El 

3. ANALYSIS 

Th e Bullet Cluster, at redshift z = 0.297 ([Tucker et al.l 
119981 ). is the product of two merging galaxy clusters. The 
less massive member has passed through the dominant 
mass with a relative velocity of > 4500 km s^^. This ex- 
treme interaction emphasizes the displacement between 
the mass of the cluster, traced by weak lensing models, 
and the hot X-ray gas, provicfijig_dramatic evidence for 
the presence of dark matter (iClowe et al. 2004). Figure[l] 
shows maps of the Bullet Cluster at optical. X-ray, and 
infrared wavelengths. The image in the left-hand panel 
combines optical maps from Magellan a nd the Hubble 
Space Telesco pe (HST) ( Clowe et^|2004D with Chandra 
X-ray data (M arkevitch et al. 200*^ The X -ray emission 
reveals hot gas that represents the majority of the bary- 
onic matter in the cluster (shown in pink). The distribu- 
tion of the dark matter is shown in blue with weak lens - 
ing contours over-plotted in green owe 
The yellow cross marks the position of the AzTEC point- 
source detection at 1.1mm, near the peak of the largest 
mass distribution. 

The right-hand panel of Figure [T] shows a 2' x 2' detail 
of the region. The color image is produced from Spitzer 
IRAC 3.6, 4.5, and 8.0 /im maps. The BLAST cen- 
troids at 250, 350, and 500 /im are shown in blue, green, 

^ Maps and PSFs for the field are available to download from 
http: / /blastexperiment.info^ 
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and red, respectively. The AzTEC 1.1 mm position is 
shown in yellow. Each cross indicates the position of 
the peak BLAST detection, with circles representing the 
FWHM of each beam. It is apparent that the relative off- 
sets between the BLAST beams is inconsistent with the 
p" pointing e r ror d et ermined in o t her B LAST06 maps 
iPascale et all (|2q08l ): iDevlin et all (|2009l ). We therefore 
consider the possibility that the centroids are shifted due 
to emission from a combination of galaxies at different 
redshifts, slightly displaced along the line-of-sight. 

3.1. Photometry 

The displacement of the BLAST centroid positions 
from the location of the millimeter source increases with 
frequency and detections are elongated in the direction 
of that displacement. These features are illustrated 
in Figure [2l showing a detailed 2'x 2' region around 
BLAST J065837-555708. The differential shift between 
bands suggests emission from other galaxies with small 
angular offsets to the line-of-sight, south of the AzTEC 
source. The central panel of Figure [2] is a composite 
of IRAC 3.6, 4.5, and 8.0 /im images. The maps reveal 
multiple cluster member galaxies that may be contribut- 
ing flux to the BLAST submillimeter measurements and 
shifting the centroid positions away from the background 
LIRG. For a galaxy contained within the Bullet Clus- 
ter (z ~ 0.3), BLAST colors sample the Rayleigh- Jeans 
side of the SED. This emission scales as S\ ^ and 
would contribute to the BLAST bands in a form con- 
sistent with the observations. The MIPS 24 /xm map is 
shown in the right-hand panel of Figure [21 with cyan cir- 
cles around > 5-ct sources. Those located within the 3-tT 
flux density contours of the BLAST emission are num- 
bered and infrared photometry extracted for each using 
the APEX pipeline is given in Table [H Though a num- 
ber of these galaxies likely contribute to the submillime- 
ter detection, sources 2, 3, and 4 clearly dominate the 
emission at 24 yum and line up with the centroid peaks of 
the submillimeter maps. Sources 2 and 3 are likel y dou- 
ble images of the lensed bac kground LIRG ( Wilso n et al.l 
I2008t iGonzalez et"al1 l2009f ) . while Source 4 is a spec- 
troscopically confirmed cluster member (|Mehlert et al.l 
1200 It) . We estimate the submillimeter flux densities of 
each of the lensed high-redshift millimeter source de- 
tected by AzTEC and the foreground dominant cluster 
member galaxy by simultaneously fitting the PSF to the 
BLAST maps at the positions of both sourc es using the 
maxim um-likelihood algorithm presented in iScott et al.l 
(|2002f ). Because the mid-infrared images of the lensed 
LIRG are closely spaced compared to the scale of the 
resolution, we adopt the millimeter position of Wilson et 
al (2008) for this component in the joint-fit. The sub- 
millimeter photometry of the LIRG and the bright clus- 
ter member galaxy are given in Table [4] and Table [5] re- 
spectively. The measured 500 //m emission is dominated 
by an association with the millimeter AzTEC source. 
However, our estimates of the shorter wavelength con- 
tributions to the flux densities at 250 and 350 /im fa- 
vor the cluster member galaxy as the dominant source 
of emission. Consequently, the observed submillimeter 
spectrum of the background LIRG is flattened across the 
BLAST bands. 

Submillimeter emission in this region is not confined 
to compact sources. IRAS maps reveal foreground cirrus 



emission of 5MJysr ^ at 100 /im. However, this struc- 
ture is smooth on scales of < 1° and therefore filtered 
out of the BLAST maps. The measured BLAST flux 
densities are adjusted to account for a small contami- 
nation by the thermal Sunyaev-Zel'dovich effect (SZE). 
This excess flux density results from the hot gas in the 
cluster boosting cosmic microwave background photons 
to higher energies through in verse Compton scattering. 
Following IWilson et all ()2008[ ). we estimate the SZE flux 
contribution to the BLAST bands at the position of the 
AzTEC detection by developing a model to describe the 
electron density distribution of the Bullet Cluster. We 
assume a two-component spherically-symmetric /3-model 
{(3 = 0.7) and a range of total cluster-mass, Mtot = 1-5- 
2.1 X 10^^ M0, distributed in a 3:1 ratio (jNussed [20M ) 
between the main and sub-clusters. Each component is 
centered at the peak of the corresponding observed X- 
ray su rface brightness distribution given bv lClowe et al.l 
(|2006f l. Provided that the electron temperature is similar 
to the X-ray gas temperature, and assuming each com- 
ponent is isothermal, we estimate electron temperatures 
^ 15.4-16.3 keV for the main cluster and ~ 6.0-7.4 keV 
for the sub-cluste r, in agreement w i th th e temperature 
map provided bv iMarkevitch et al.l (|2002l ). Relativistic 
corrections w ere applied using the analytical approxima- 
tion given bv lltoh et al.l (| 19981 ). The results of these sim- 
ulations and our revised source flux densities are given 
in Table m We have also used this model to estimate the 
SZE at the position of the bright cluster member galaxy 
and these results are shown in Table [5l 

3.1.1. Photometric Redstiifts 

The shape of the spectral energy distribution (SED) 
in the submillimeter regime can be used to estimate the 
redshift of the source. The BLAST bands were chosen 
to optimize constraints on the FIR peak of the bright 
submillimeter galaxy population at 1 < z < 4. A Monte 
Carlo photometric redshift method has be en developed 
(jHughes et al.l 120021: lAretxaga et"aLl I2003D to take into 
account constraining prior information such as a range 
of local SED templates and the evolving luminosity func- 
tion to z « 2. 

In accordance with this method, we generate SEDs 
for a catalog of galaxies based on an evolutionary model 
for the 60 /im luminosity function that flts the observed 
850 /im number counts (e.g. luminosity evolution oc (1 -I- 
z)^ for z < 2, and no evolution at z > 2) and covers 
a simulated area of lOdeg^ Template SEDs are drawn 
at random from a library of local starbursts, ULIRGs, 
and AGN to provide FIR-radio fluxes. The SEDs cover 
a wide range of FIR luminosities (10^"-10^^ '^ Lq) and 
temperatures (25-65 K). The photometric and calibra- 
tion uncertainties of the BLAST and AzTEC measure- 
ments are included in the fluxes of the mock galaxies. 
We adopt a conservat ive estimate for the amplification of 
the s ource {A = 30) (jWilson et al.l 120081: [Gonzalez et all 
[2009"), and reject mock galaxies with fluxes below the de- 
tection threshold of the survey. The redshift probability 
distribution of the source is then calculated as the nor- 
malized distribution of the redshifts of the mock galaxies 
in the catalog, weighted by the likelihood of identifying 
the colors and fluxes of each mock galaxy with those of 
the source. 
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Figure [3] shows the results of this analysis indicat- 
ing a photometric redshift Zphot = S.Qto'J- This esti- 
mate is consistent with the IR-color redshift estimation 
^ 2.7 ± 0.2) for the most likely IRAC counter- 
part (jWilson et al.ll2008l : iGonzalez et al.ll2009D . The ac- 
curacy of the technique has been experimentally tested 
for galaxies detected in at least 3 bands in the rest- 
frame radio to FIR interval, yielding a typical l-cr error 
((-grtiot--gBpec)^)^/^ ~ 0.25 over the 0.5 < z < 4.0 regime 
(jAretxaga et al.ll2007| ). 



3.1.2. Intrinsic Properties 

With an estimate of its redshift, we can derive prop- 
erties intrinsic to the submillimeter galaxy. Figure |4] 
shows a single-temperature modified blackbody SED of 
the form 

= N{v/v^fB,{T), (1) 

where is the flux density, N is the amplitude, [3 is 
the dust emissivity index, is fixed at c/250/im, and 
B^(T) is the Planck blackbody radiation function for a 
source with temperature T. The function is fit to the 
data using a minimization routine, after subtracting 
the contribution from SZE and an estimate of the contri- 
bution from the brightest mid-infrared cluster member 
galaxy, as discussed previously. The grey curves indicate 
the l-tJ range of temperature and amplitudes for the fit 
estimated from 200 Monte Carlo simulations. The value 
of /3 is fixed here at (3 = 1.5, but the fits give very similar 
values of T and Lfir if we use (3 = 2.0. The best-fit 
temperature of the distant lensed source is 32 K. This 
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). This spectral information also enables a color 



correction to the flux densities which acco unts for the 
profil es and widths of the BLAST filters (|Truch et all 
l2008f ). These color corrected flux densities are given 
in Table ID Lpj^ is the FIR luminos ity, calculated i n 
the wavelength range A = 8-1000 ^m (iK ennicufrBllQQSf ). 
for the observed flux densities without considering grav- 
itational lensing effects. Using the previously assumed 
lower limit for the amplitude of the lensing magnifica- 
tion {A > 30), we derive an upper limit for the intrin- 
sic FIR luminosity of Lp^^ < 10^'^ Lq. This value is 
consistent with esti mates presented by the AzTEC team 
(jWilson et al.ll2008l) as well as from extrapolation of the 
24 /xm luminosity (jGonzalez et al.l [2009) , and qualifies 
the source as a LIRG. This intrinsic luminosity implies 
an upper limit on the star formation rate of 200 Mq yr~^ 
(|Kenmcutli[l99l . The source is relatively faint com- 
pared to typical submillimet er-selecte d galaxies wi th sim- 
ilar dust temperature s ( B lain et al.l [2003: Kova cs et al.l 
[2001 IDve etaHllffli ). and is more consistent with the 
luminosity-temper ature rela tionship obs erved in IRAS- 
selected samples (iSoifer et a l. 1991; Bl ain et all 120031 : 
iChapman et al.ll2003t iChapin et al.ii2009f) . 

4. CONCLUSIONS 

BLAST submillimeter maps of the Bullet Cluster have 
revealed bright e mission near the peak of the weak 
lensing contours (' Clowe et al.l I2006D . The position of 
the detection coincides with the AzTEC 1.1mm source 



MMJ065837-5557.0 and a multipl y-imaged LIRG de- 
tected in the Spitze r IRAC maps (jWilson et al.l 120081 : 
IGonzalez et al.l [20090 . The systematic offsets in peak 
emission between the BLAST bands motivates the likely 
interpretation that multiple sources are confused at 
BLAST resolution. We have performed a simultane- 
ous fit to the brightest two infrared galaxies in the field 
in order to estimate the properties of the background 
LIRG. The resulting color corrected flux densities for 
the the source are S'250 = 94 ± 30, S'350 = 96 ± 27, and 
5500 = 110±21mJy. 

The BLAST colors sample the peak of the source spec- 
trum, and 250-1 100 /im flux densities place strong con- 
straints on its photometric redshif t. Using a Monte 
Carlo photometric red shift method ([Hughes et al.ll2003 
lAretxaga et al.l l2003f l. we have calculated a redshift 
-Zphot — 2.9lg'3. '^^^ redshift probability distribution 
is inconsistent with the redshift of the Bullet Cluster 
(z « 0.3), ruHng out the possibility that the miUimeter 
AzTEC source is associated with the foreground mass. 
Assuming the best-fit photometric redshift we have de- 
rived the observed FIR luminosity of the source, in the 
absence of lensing; to be = 2 x 10^^ L©. The 

proxim ity of the source to a critical line in the lensing 
model ([Bradac et al.l 120061 : IGonzalez et al.ll2009f l puts a 
lower limit A > 30 on the amplitude of magnification 
due to gravitational lensing by the cluster. This imposes 
an upper limit on the intrinsic luminosity of the SMC, 
Lpjp < 10^^ LiT). This luminosit y is consistent with 
estima tes from iWilson et all (|2008f ) and IGonzalez et all 
(|2009( ). and characterizes the source as a LIRG, a mem- 
ber of the fainter population of SMGs that dominate the 
number counts but are difficult to detect with the limited 
resolution of existing submillimeter telescopes. 

The challenges in interpreting the BLAST maps 
demonstrate the necessity of including complementary, 
higher resolution data sets in the analysis. BLAST is a 
precursor to the Spectral and Photometric Imaging Re- 
ceiver (SPIRE ) on the recently launched Herschel Space 
Observatory ( Griffin et al.ll2006[ ). With nearly twice the 
resolution, and gains in sensitivity, SPIRE will be able to 
provide significantly better centroiding accuracy. Never- 
theless, access to infrared and millimeter data sets will 
still be integral to a more complete understanding of 
these crowded systems. 
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ing Research Council (NSERC) of Canada, and the UK 
Science and Technology Facilities Council (STFC). We 
thank Tony Mroczkowski for his help with the SZE sim- 
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West Grid computing resources. This research also made 
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Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA, and SAOImage DS9, 
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Fig. 1. — Left panel: T he background image combines optical maps from Magellan and HST HClowe et al.ll200^ with Chandra X-ray data 
(jMarkcvitch ct al."2003). The X-ray map reveals the hot gas which represents the majority of the baryonic matter in the cluster (shown in 
pin k). The distribution of the dark matter is shown in blue, with weak lensing contours over-plotted in green. Image based on Figure 1 in 
Oowe et al.i (,2006 ). The yellow cross indicates the position of the AzTEC detection. Right panel: A 2'x 2' detail based on IRAC 3.6, 4.5, 
and 8.0 fim maps of the region. The BLAST 250, 350, and 500 /im positions are shown in blue, green, and red, respectively. The AzTEC 
1.1 mm position is shown in yellow. The crosses indicate the centroids of the detections and the circles re present the FWHM of each beam. 
The white contour denotes the critical curve corresponding to infinite magnification based on HST data IBradac et al.ll2006l) . 
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TABLE 1 

Millimeter properties of AzTEC source 
MMJ065837-5557.0. 



h'lux Density 

RA Centroid Position 

Dec Centroid Position 



13.5 ± l.OmJy 

ei^ss^sTfai {J2000) 

-55°57'1"5 (J2000) 



Note. — 1.1mm flux density and cen- 
troid positio n of MMJ065837-5557.0 from 



IWilson et all | [200l) . 
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Fig. 2.— Detailed 2'x 2' region around BLAST J065837-555708. 
Tlie BLAST 250, 350, and 500 fim positions are shown in blue, 
green, and red, respectively. Crosses indicate the centroids of the 
measurements and corresponding solid circles represent the FWHM 
of the beam. Dashed lines show 3-cr contours for each BLAST de- 
tection. Left-hand panel: HST 606 nm map. Small green circles 
indicate the position s of spectroscopically verifled cluster members 
IIMehlert et al.ir20011) . Many of these foreground galaxies may be 
contributing flux to the submillimeter measurement and shifting 
the centroid positions away from the position of the background 
LIRG. Middle panel: Red, green, and blue composite image made 
from IRAC 8.0, 4.5, and 3.6 fim maps, respectively. Again, cluster 
members are indicated with small green circles. 24 fim counter- 
parts detected with a signiflcance > 5-cr are marked in cyan. The 
infrared maps reveal a number of bright galaxies near the peak of 
the 350 fj,m emission. Right-hand panel: MIPS 24 ^m map. Cyan 
circles indicate > 5-(t sources. Sources within the 3-cr BLAST con- 
tours are numbered and the extracted infrared photometry is given 
in Table |3] Sources 2 and 3 are likely split images of the lensed 
background LIRG associated with the 1.1mm detection. They co- 
incide with the peak of the 500 /im emission, although the shorter 
wavelength peaks fall closer to source 4. This lower redshift source 
is likely to be brighter at short BLAST wavelengths, contributing 
signiflcantly to the emission in those bands. Sources 1, 5, and 6, 
which also appear at 24 /im may additionally contribute at BLAST 
wavelengths. The elongation of the contours could be explained by 
confusion of these sources in the lower resolution BLAST maps. 



1.000 



0.100 



0.010 



0.001 



■^250 = 94 +/- 30 mjy 

Sjso = 96 +/- 27 mJy 

S500 =110 +/- 21 mJy 

S,,oi, = 13.5 +/- 1.0 mJy 
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I . 



32K 
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Lj;i < 1x10'' Le 
SFR < 200 Mg yr- 



BLAST_J065837 



1000 



100 

Wavelength [fj,m] 



Fig. 4. — Spectral Energy Distribution of 
BLAST J065837-555708. The black curve is the best-fit 
model to the BLAST and AzTEC measurements, shown on the 
plot with 1-(T error bars. These flux densities have been corrected 
for the contribution from SZE. The grey curves indicate the l-cr 
range of temperature and amplitudes for the flt. The spectral 
index is fixed at /3 = 1.5. The spectral information enables a color 
correction to the measured flux densities which accounts for the 
proflles of the BLAST fllters. These color corrected values are 
displayed on the plot, ipjf^ represents the far-infrared luminosity 
based on the fit to observations, while iJJ?j'j^ is an estimate of 
the intrinsic far-infrared luminosity of the source, assuming an 
amplification of A > 30 due to gravitational lensing. An upper 
limit on the star formation rate is derived based on this luminosity 
HKennicuttllTQOai) . 



o 
n 
o 




Fig. 3. — The redshift probability distribution based on 250- 
1100 /.tm photometry at the millimeter centroid position. The un- 
certainty in lensing magnification results in a longer high-redshift 
tail. 
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TABLE 2 

Peak Positions of BLAST J065837-555708. 



A 

[^tm] 


Before Correction 
Centroid RA Centroid Dec 
[h:m:s] [d:m:s] 


After Correction 
Centroid RA Centroid Dec 
[h:m:s] [d:m:s] 


250 


06:58:35.3 


-55:57:08 


06:58:36.7 


-55:57:33 


350 


06:58:36.2 


-55:56:58 


06:58:36.6 


-55:57:22 


500 


06:58:36.0 


-55:56:44 


06:58:37.3 


-55:57:09 



Note. — Positions of the peaks in each band associated with 
BLAST J065837-555708. The positions are given before and after the 
correctional offset applied based on stacking the 250 fim map against 
24 /xm sources. 



TABLE 3 

Infrared photometry for 2AfiM sources associated with BLAST J065837— 555708. 



ID 24 i^m R A 


24 ^m Dec 


S3. 6 


S4.5 


S5.8 


Sg.o 


S24 


[Hours] 


[Degrees] 




[MJy] 


[MJy] 


[MJy] 


[MJy] 


1 6:58:39.69 


-55:56:56 


18.6 ± 0.2 


26.2 ± 0.3 


70.6 ± 1.1 


40.0 ± 1.9 


252.7 ± 37.2 


2 6:58:37.99 


-55:57:01 


19.1 ± 0.2 


25.7 ± 0.3 


76.5 ± 1.1 


60.5 ± 1.9 


596.2 ± 34.9 


3 6:58:37.21 


-55:57:06 


28.5 ± 0.2 


49.0 ± 0.3 


69.2 ± 1.1 


86.9 ± 2.0 


840.4 ± 36.4 


4 6:58:37.40 


-55:57:24 


309.6 ± 0.5 


269. 5± 0.5 


190.4 ± 1.3 


528.3 ± 2.1 


1482.0± 34.8 


5 6:58:35.62 


-55:57:55 


72.0 ± 0.2 


47.4 ± 0.3 


60.2 ± 1.1 


25.6 ± 1.9 


306.4 ± 37.4 


6 6:58:34.09 


-55:57:54 


32.5 ± 0.2 


23.0 ± 0.3 


41.9 ± 1.1 


18.4 ± 1.9 


285.0 ± 36.5 


Note. — IR photometry corresponding to the 24 rim sources located within the BLAST 3-cr contours 


which have a significance of > 5- 


















TABLE 4 










SUBMILLIMETER FlUX DENSITIES OF MM,I065837- 


'5557.0. 




A 


Measured Flux 


SZE 


SZE Corrected Flux Color Corrected Flux 


[^tm] 


[mJy] 


[mJy] 


[m 


■Jy] 


[mJy] 




250 


91 ± 27 


0.15 ±0.03 


91 ±27 


94 ±30 




350 


95 ± 24 


1.1 ±0.2 


94 ±24 


96 ±27 




500 


121 ± 16 


10.6 ±1.7 


110 ±16 


110 ±21 




Note. 


— Flux densities in each band associated with BLAST J065837— 555708 at 



the position of the AzTEC detection. The contribution from SZE is calculated at 
the position of the AzTEC point source, based on the theoretical model outlined in 
[Wilson et al. ( 2008). This small contamination is subtracted from the measured fluxes 
for further analysis. The monochromatic flux densities based on the best-fit SED are 
given in the final column. 



TABLE 5 

SUBMILLIMETER FlUX DENSITIES OF THE DOMINANT CLUSTER MEMBER GALAXY. 



A 

[mH 


Measured Flux 
[mJy] 


SZE 

[mJy] 


SZE Corrected Flux 
[mJy] 


Color Corrected Flux 
[mJy] 


250 


129 ± 26 


0.13 ±0.03 


129 ±26 


131 ±30 


350 


104 ± 23 


1.0 ±0.2 


103 ±23 


106 ±27 


500 


49 ± 18 


9.4 ±1.5 


40 ±18 


38 ±18 



Note. — Submillimeter flux densities in each band associated with the dominant 
24fim cluster member galaxy at RA (J2000) = 6:58:37.40, Dec (J2000) = -55:57:24, 
identified as Source 4 in Figure [2] and Table [3] The columns are analogous to those in 
Tabled 
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